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lO ■ One method to search for particle dark matter is to hunt down its annihilation products. 

' In the Solar System, three potential types of signals of annihilation have been identified: 
neutrinos and 7-rays from the Sun, and neutrinos from the Earth. Each of these signals 

1/^ ' depends sensitively on the orbital evolution of dark matter once it becomes bound to the Solar 

' System. I will review progress on characterizing these signals based on recent improvements 

0^ . in the determination of the properties of the bound dark matter population. 
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■ 1 Introduction 

■ Many lines of evidence point to the existence of dark matter in the universe, but its composition 
remains a mystery^H. Although other particle physics (e.g., the axion^^ and astrophysical (e.g., 
primordial black holes^ candidates exist, by far the most popular candidate for dark matter is 
one or more species of weakly interacting massive particle (WIMP). Many proposed extensions 
to the Standard Model naturally produce such particles^l. 

Numerous experimental methods are being used to determine the nature of dark matter. 
Direct detection experiments are designed to measure the energy of recoiling nuclei after they 
have been hit by astrophysical WIMPsEI. Byproducts from WIMP annihilation in dark matter 
concentrations throughout the Galaxy may be probed with cosmic ray experiments ^. 

The two brightest sources of annihilation products in the Solar System should be the SuiP 
and EarthEl. WIMPs generically have small interaction cross sections with nuclei, and as such 
can be captured and settle into dense cores in the potential wells of celestial bodies. Thus, 
WIMP densities n in these bodies can be many orders of magnitude above the local Galactic 
density, which is important since the annihilation rate per unit volume dT/dV oc n^. WIMP 
annihilation in the Earth may be detectable since it is a nearby source, while the Sun's far 



deeper potential well (and a correspondingly higher WIMP capture rate) compensates for its 
much greater distance from terrestrial detectors. 

There are three possible signatures of WIMP annihilation in these bodies. The only an- 
nihilation products which may escape from the interiors of the Sun and Earth are neutrinos. 
The muon neutrinos may be observed in terrestrial neutrino telescopes (e.g., IceCubeEl, Antares 
1^ via the Cerenkov radiation of neutrino-induced muons in and around the detector volume. 



These signals may be distinguished from backgrounds by their directionality and energy spec- 
trum. Strauszim suggested that WIMPs annihilating just outside of the Sun may be visible in 
7-rays ("near-solar 7-rays"). Even though the annihilation rate is much lower outside the Sun 
than in its core, Strausz suggested that since the Sun is expected to emit few photons at high 
energies during its evolution, the backgrounds to the near-solar 7-ray signal would not be severe. 

In the next sections, I will describe the standard assumptions used in the calculations of 
these signals, and show how the calculations change when the details of the evolution of WIMPs 
bound to the Solar System are included. 

2 Neutrinos from the Sun 

The upper limit on the neutrino flux from WIMP annihilation in the Sun^SHSl currently places 
the tightest constraint on the WIMP-proton spin-dependent elastic scattering cross section <Tp^. 
The IceCube e2cperiment is expected to have at least an order of magnitude better sensitivity 
once complete^ as will the proposed KMSNeT^l^. These experiments are projected to probe 
significant parts of the beyond Standard Model phase space, although this statement depends 
critically on the assumptions used to calculate event rates. 

Schematically, the event rate of neutrino-induced muons in neutrino telescopes can be de- 
scribed by 

A'^^ oc r X (neutrino physics) x (detector properties) , (1) 

where T is the total annihilation rate of WIMPs in the Sun and "neutrino physics" encompasses 
details on the annihilation branching fractions, propagation through the Sun and to the tele- 
scope, and interactions near the detector. The standard calculation of the annihilation rate 
includes the following assumptions: (1) The thermalization is nearly instantaneous. (2) All 
WIMPs that are captured in the Sun thermalize. (3) The local Galactic WIMP distribution 
function is spatially smooth, and with velocities distributed as Maxwellian in an inertial Galac- 
tocentric frame (and with a one-dimensional velocity dispersion a = 155 km s~^). I will touch 
on point (3) in Sec. 13.11 but describe the problems with the first two assumptions here. 

Assumption (1) is violated due to the finite optical depth in the Sun to WIMPs, r, and 
the finite energy a WIMP may lose to a solar nucleus in each encounter, Q. Since the Sun 
must be optically thin to WIMPs (deduced from existing constraints on the elastic scattering 
cross sections), WIMPs may make many passages through the Sun before rescattering, with the 
typical time between scatterings scaling as tr ~ -fx/''"' '^here is the WIMP orbital period. In 
addition, it generically takes many scatters for a WIMP to thermalize because 

Q ~ mAvl^^, (2) 

but the WIMP orbital energy is 

E mxt;2^^(i?0/a), (3) 

where uia is the mass of a solar nucleus with atomic number A, is the WIMP mass, Vesc 
is the escape speed from the surface of the Sun, and a is the semimajor axis. A thermalized 
WIMP should have an energy 

Etherm <■ -m^vl^c- (4) 
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Figure 1: Suppression of the annihilation rate due to the gravitational interactions of WIMPs with Jupiter in 
the cases that (left) spin-independent and (right) spin-dependent elastic scattering dominate in the Sun. Fa is 
the annihilation rate estimated from my simulations, and is the annihilation rate calculated using standard 

assumptions. Figure from Peter 



As the WIMP mass increases, more scatters are required for the WIMP to sink to the center 
of the Sun. Furthermore, as the WIMP mass increases, the median semimajor axis of the orbit 
to which a Galactic WIMP is captured increases, which imphes that heavy WIMPs must lose 
more specific energy than lighter WIMPs to thermalize. Thus, the median thermalization time 
for the solar captured WIMPs increases as the WIMP mass increases. If spin-dependent (spin- 
independent, denoted SI) scattering dominates in the Sun, the median captured WIMP will 
require the age of the Solar System to thermalize if cXp^ ~ 10~^^ cm^ (ap^ ~ 10^^^ cm^) if 
= 100 GeV, but « IQ-"^^ cm^ {a^^ lO"'^^ cm^) if = 10 TeV. This is described in 
more detail in PeterlJ^. 

Assumption (2) is violated due to the presence of planets in the Solar System, whose grav- 
itational torques can affect orbits in two generic ways. They can eject the WIMPs from the 
Solar System, meaning that those WIMPs will never thermalize; or they can alter the perihelia 
of the orbits so that WIMPs either pass less frequently through the Sun or only in the outskirts 
in the Sun where the optical depth is much lower. In the latter case, thermalization may be 
significantly delayed, even beyond t©. 

I used a set of WIMP orbit simulations in a simplified solar system consisting of only the Sun 
and Jupiter (originally performed to determine the phase space density at the Earth of WIMPs 
bound to the solar system^l^l) iq determine which effects dominated for a given initial semimajor 
axis. I found that WIMPs initially scattered onto orbits with a > 2.6 AU were mostly ejected 
from the Solar System unless the thermalization time tt was less than the time required for the 
gravitational torque from Jupiter to lift the WIMP perihelion from the Sun (~ 1000 years). 
WIMPs with 1.5 AU < a < 2.6 AU had rescattering times tr ~ 300P-^/t. This is longer than if 
the Sun were in isolation because secular and mean-motion resonances pull the WIMP perihelia 
out of the Sun for extended periods of time. WIMPs with initial a < 1.5 AU had thermalization 
times largely unaffected by planetary torques. 

The effects on the annihilation rate of violating assumptions (1) and (2) are shown in Fig. [U 
in which I plot the ratio of the annihilation rate estimated from my simulations to the standard 
calculation as a function of WIMP mass. Lines for several values of the elastic scattering cross 
section are show. I fix the annihilation cross section to (av) = 3 x 10^^^ cm^ s~^, and use 



a conventional halo model for the local dark matter distribution. In the left (right) panel, I 
show the suppression if spin-independent (spin-dependent) scattering dominates in the Sun. For 
cTp^ > lO"^'^ cm^ ('^p^ > 10~^^ cm^), there is no suppression because even the orbits with 
a > 2.6 AU thermalize before gravitational torques from Jupiter can significantly affect the 
orbits. WIMPs with a < 2.6 AU also thermalize rapidly. For the uppermost line in each panel 
of Fig. 1, the annihilation rate is suppressed due to the ejection of WIMPs with a > 2.6 AU. 
However, WIMPs with a < 2.6 AU thermalize on timescales less than t^. The suppression is 
greater at higher WIMP masses because an increasing fraction of WIMPs are initially captured 
onto long orbits, and is greater if spin-dependent scattering dominates in the Sun because the 
target nucleus on which the WIMPs scatter is much lighter than if spin-independent scattering 
dominates. This reduces the typical energy loss per scatter (Eq. ([1])). The middle line in each 
panel shows suppression due to ejection, but at this point the capture and annihilation rates 
drop out of equilibrium in the Sun. Whereas P oc C, where C is the capture rate of WIMPs in 
the Sun, if the rates are in equilibrium (with a coefficient that depends on whether WIMPs are 
self-annihilating or not), F oc if the rates are far out of equilibrium. In the lowest lines, the 
cross section is small enough that the thermalization time of the 1.5 AU < a < 2.6 AU WIMPs 
is longer than t©. For cross sections much lower than these values, the thermalization time for 
WIMPs with a < 1.5 AU exceeds Iq, causing a near-total suppression of the annihilation rate. 

In summary, the next generation of neutrino telescopes will have far lower sensitivity to high 
mass {m^ > 1 TeV) WIMPs than the standard annihilation calculation would suggest. 

3 Neutrinos from the Earth 

The center of the Earth is an attractive target for WIMP searches due to its proximity to 
neutrino telescopes. However, there are two complications in searching for annihilation in the 
Earth. First, direct detection experiments constrain the spin-independent WIMP-proton cross 
section much more tightly than the spin-dependent cross section. Since the Earth only has trace 
elements that may interact via spin-dependent channels with WIMPs (the dominant isotopes 
in the Earth, ^^Fe, ^^O, and ^^Si, may only have spin-independent WIMP interactions), the 
capture rate is far more constrained than the capture rate of WIMPs in the Sun (in which the 
dominant species, hydrogen, may have spin-dependent interactions). 

Second, the Earth's potential well is shallow — the escape speed from the center of the Earth 
is Vesc ~ 15 km s~^ compared to the Sun's ~ 1000 km s~^. Typical halo WIMPs have speeds an 
order of magnitude higher than the escape speed from the Earth, which makes them kinemat- 
ically difficult or impossible to capture in the Earth unless the WIMP mass is nearly the mass 
of one of the nuclear isotopes in the Earth^^Sl^ jf j^^^ > 4Q0 GeV, halo WIMPs are impossible to 
capture; the only WIMPs the Earth may capture are those already bound to the Solar System. 

Thus, in order to predict event rates at neutrino telescopes (or to derive a WIMP annihilation 
rate from a signal), the Solar System's bound WIMP population must be characterized. Using 
detailed balance arguments based on the gravitational scattering of WIMPs by planets, Gould 
claimed that the annihilation rate could be accurately derived using the phase space density 
of WIMPs far outside the gravitational sphere of influence of the Solar System (the "free space" 
density). More recently, Lundberg & Edsjo^^ solved a diffusion equation for WIMPs in the 
Solar System, taking into account the loss of WIMPs in the Sun due to thermalization. They 
found that the phase space density of bound WIMPs at the Earth was greatly reduced if the 
Sun were infinitely optically thick to WIMPs than if it were described as a gravitational point 
source. Damour & Krauss^ described a population of WIMPs that could be captured by elastic 
scattering in the Sun, and survive for long times in the Solar System due to a secular resonance 
that forced the perihelion distance to the center of the Sun to oscillate. If this population survived 
the lifetime of the Solar System, it could enhance the annihilation rate in the Earth by almost two 
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Figure 2: Neutrino-induced muon fluxes from WIMP annihilation in the Earth above an energy threshold Ej^ = 1 
GeV for water/ice neutrino telescopes for (left) unbound halo WIMPs and (right) the maximum expected from 
both bound and unbound dark matter, in scans of the MSSM. The solid line represents an optimistic detection 
threshold for the IceCube telescope The open circles represent neutralino models with spin-independent cross 
sections above the 2006 CDMS exclusion curve , the triangles represent models with cross sections between the 
2006 and the present limits, and models marked with blue squares are consistent with current bounds on the 

elastic scattering cross sections Figures from Peter 



orders of mae nitude if 60 GeV < < 130 GeV^. Annihilation rates are usually calculated 
using the free space distribution function, the purely unbound WIMP distribution function 
(which may be calculated at the Earth using Liouville's theorem and Galilean transformations), 
or Lundberg & Edsjo's phase space density. 

Dynamics in the Solar System are far more complicated than can be described using the 
diffusion or semi-analytic treatments used in the work cited above to estimate the bound WIMP 
population. The only way to truly determine the bound WIMP density in the Solar System 
is to simulate a statistically significant number of orbits. As a first step, I simulated bound 
WIMP orbits in a solar system consisting only of the Sun and Jupiter. I simulated both orbits 
bound to the solar system by elastic scattering in the Sunlit and WIMPs that are captured by 
gravitational interactions with JupiterEI. I found that the WIMPs bound by the first method 
have shorter lifetimes than assumed by Damour & Krauss, and thus have a significantly lower 
density. WIMPs that are gravitationally captured to the Solar System (the population explored 
by Gould and Lundberg &: Edsjo) have a slightly lower density than found by Lundberg & Edsjo. 

To estimate the neutrino-induced muon flux at the surface of the Earth, I assumed that the 
WIMP was the neutralino in the minimal supersymmetric standard model (MSSM), and used 
the DarkSUSY software package^^^l to make scans of MSSM parameter space. The results are 
shown in Fig. [2j Fluxes due to unbound Galactic WIMPs only are shown on the left; those 
due to unbound WIMPs and the maximum possible bound WIMP density are on the right. 
The solid line indicates an optimistic detection threshold for the IceCube neutrino telescope. 
The squares mark models that are still allowed by direct detection experiments. I find that the 
bound WIMPs only marginally increase the fiuxes at the Earth, and that allowed models lie 
well below threshold. 

There are two caveats to this pessimistic conclusion. First, I have only shown the fiux of 
through-going muons; since IceCube is large compared to typical muon path lengths through 
ice, there could be a large number of contained events. However, this is likely to increase the 
event rate above the through-going event rate only by factors of several to about a factor of ten. 



This is still insufficient to boost the muon ffux above threshold. 

Second, the regions of WIMP phase space that are most accessible to capture are not ac- 
cessible in the simplified solar system. This part of phase space can only be populated by close 
encounters between WIMPs and the inner planets. Therefore, more sophisticated simulations 
are required to accurately predict the signal of WIMP annihilation in the Earth. 

3.1 Everything Is Not Lost: Dark Matter Disks and Detection 

Although the prospects for detecting WIMP annihilation in the Earth are grim, these conclu- 
sions are drawn based on fairly strict assumptions on how dark matter is distributed in the 
Galaxy, and what its phase space structure is. It is almost always assumed that the dark matter 
is distributed almost smoothly in a nearly spherically symmetric Navarro-Frenk- White haloes 
about the Galactic center. Moreover, it is assumed that the halo is at rest in an inertial Galacto- 
centric frame, and that the velocity distribution can be described locally by a Maxwellian with a 
velocity dispersion a = vq/\/2, where the stellar disk has a rotation speed of Vq = 220 km s~^. 
These assumptions are based on N-body simulations of dark matter in roughly Milky Way-sized 
halos. 

However, the Sun lies in part of the Galaxy that is dominated by luminous matter, not dark 
matter. Any local dark matter parameters that are derived from simulations should be derived 
from simulations that include luminous matter, too. Recently, Read et al. have analyzed 

simulations of disk galaxies in Milky Way- mass halos, in both idealized and cosmological con- 
texts. They find that, in general, massive satellite galaxies are dragged into the spiral galaxy 
disk plane and dissolved, yielding a "dark disk" of dark matter. This component is in addition 
to the conventional, more spherical dark matter halo. The properties of the dark disk depend 
sensitively on the accretion history of the galaxy, but the "median" galaxy in their samples had 
a dark disk with a rotation velocity lag with respect to the Sun of 50 km s^^, and a velocity 
dispersion of <t ~ 50 km s~^, and a mass density in the mid-plane similar to that of the halo. 

This dark disk is important for WIMP detection because the typical WIMP speed with 
respect to the Solar System is much lower than for WIMPs in the halo. The capture probability 
in the Solar System increases dramatically with decreasing relative speed. In work lead by 

inol 

Tobias Bruchi^, we found the consequence was that the neutrino flux from WIMP annihilation 
in the Earth could be increased by two or three orders of magnitude assuming the median dark 
disk, and the neutrino flux from the Sun could be boosted approximately an order of magnitude. 
Even for the least signiflcant dark disk found in the Read et al. samples, the neutrino flux from 
the Earth and the Sun could be enhanced by a factor of order unity. Therefore, the dark disk 
may boost the neutrino flux for WIMP models consistent with experiments (farther) above the 
detection threshold for IceCube. 

There are two major sources of uncertainty in the dark disk enhancement of the annihilation 
rate. First, the properties of the dark disk are highly uncertain. The dark disk is likely to be 
difficult to probe with stellar dynamics, and any stellar debris from satellite destruction in the 
disk will also be difficult to disentangle from other stellar populations. Dark matter experiments 
are the most sensitive probes of the dark disk. The energy spectrum of events in direct detection 
experiments should yield some constraints on dark disk properties'^, but such constraints will 
take years (and an actual detection of dark matter in more than one experiment!) to materialize. 

Second, whether a high mass WIMP is detectable or not in neutrinos from the Earth depends 
sensitively on the phase space density of WIMPs bound to the Solar System. This is illustrated 
in Fig. 3 of Bruch et al^^, which shows that for, e.g., a 400 GeV WIMP, the annihilation rate is 
about three orders of magnitude higher for Gould's free space phase density than for the phase 
densities found in my simulations. Therefore, more sophisticated versions of my simulations, 
including more planets and more realistic planet orbits, will be required to determine the bound 



dark matter phase space density to sufficient precision to estimate WIMP parameters from 
neutrino telescope data. 



4 Near-Solar 7-Rays 

Strausz^H suggested that WIMP annihilations occuring just outside the Sun may be visible in 
7-rays. Although the dark matter density is far lower outside the Sun than at its center, it is 
higher right outside the Sun than far outside the Sun's gravitational sphere of influence. This is 
due to several effects. First, gravitational focusing increases the unbound WIMP density deep in 
the Sun's potential well. Second, in the absence of kinematic suppression, WIMPs are captured 
onto bound orbits with a semimajor axis distribution d log A^/d log a = —1. Third, massive 
WIMPs (much heavier than the typical solar nucleus) lose energy during each encounter with a 
solar nucleus. WIMPs initially captured onto barely bound orbits will spend time just outside 
the Sun as they thermalize. Moreover, the Sun does not produce 7-rays at its surface as a result 
of stellar evolution. Strausz suggested the Milagro, an air shower array used to detect cosmic 
rays, should be able to detect these near-solar 7-rays. However, it was only able to place an 

on 

upper limit to the near-solar 7-ray flux . 

Recent work suggests that the WIMP-generated 7-ray signal should be virtually unde- 
tectable. Hooper 1311 and Sivertsson & Edsjo^ have redone Strausz's calculations, and find 
that the 7-ray flux should be many orders of magnitude lower than Strausz's predictions. Using 
Monte Carlo realizations of the thermalization process, Sivertsson & Edsjo predict a flux of 
~ 10"'^ km~^ yr~^ if 0"^^ = 10~^^ cm^ and = 1 TeV. Using my own simulations, I have 
also estimated the expected 7-ray flux, and find similar results, although the flux is suppressed 
even further for high mass WIMPs due to the gravitational effects described in Sec. [2j It is not 
clear what the error in Strausz's calculation is, but current calculations suggest that the flux of 
near-solar 7-rays at the Earth should be tiny. 

In addition, the Sun is more luminous 7-rays than initially postulated by Strausz. Cosmic 
rays protons can produce pion showers in the chromosphere of the Sun^SH]^ which produce 7-rays 
by TT^ — > 77. Seckel et al. estimated that the total 7-ray flux from this process should be 
~ 10^° km~^ yv~^ (~ 10~^ cm~^ s~^) above 1 MeV, and about an order of magnitude less above 
1 GeV. 

Inverse-Compton scattering of solar photons on cosmic ray electrons can also produce a 
significant halo of 7-rays around the Sun. Orlando & Strong^' find a detection of this signal 
in EGRET data, a flux of ~ 3 x lO"'^ cm'^ s"^ in the 100 - 300 MeV energy band, which is 
expected to be stronger than the cosmic ray proton-induced 7-ray flux by factors of several. The 
Fermi Gamma-ray Space Telescope (Fermi) should be exquisitely sensitive to these cosmic ray- 
induced 7-rays from the Sun, and the signal should completely drown out any WIMP-induced 
signal from the Sun. 



5 Conclusion 

In this paper, I have described recent updates to the estimates of the solar and terrestrial neutrino 
and near-solar 7-ray fluxes from WIMP annihilation at the centers of the Sun and Earth due 
to new estimates of the orbital evolution of WIMPs in the Solar System. Some of the effects I 
described may change the event rates by orders of magnitude. I would like to emphasize that in 
order to understand the particle physics of WIMPs using new astrophysical data sets (whether 
it be neutrinos or 7-rays from Solar System sources or 7-rays from the Galactic center), it is 
necessary to understand the astrophysical properties of WIMPs: their distribution throughout 
the Galaxy, on both Solar System and ~ kiloparsec scales; and the velocity structure. While 



there has been a lot of progress in constraining the astrophysical properties of WIMPs, there 
are still many large uncertainties in these properties. 
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